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ABSTRACT 

We present sub-arcsecond observations toward the massive star forming region G75. 78 -1-0.34. We used the 
Very Large Array to study the centimeter continuum and H2O and CH3OH maser emission, and the Owens 
Valley Radio Observatory and Submillimeter Array to study the millimeter continuum and recombination lines 
(H40a and H30a'). We found radio continuum emission at all wavelengths, coming from three components: 
(1) a cometary ultracompact (UC) H II region with an electron density ~3.7 x 10"* cm"-*, excited by a BO type 
star, and with no associated dust emission; (2) an almost unresolved UCH II region (EAST), located ~6" to the 
east of the cometary UCH II region, with an electron density ~1.3 x lO*" cm""*, and associated with a compact 
dust clump detected at millimeter and mid-infrared wavelengths; and (3) a compact source (CORE), located 
~2" to the southwest of the cometary arc, with a flux density increasing with frequency, and embedded in a 
dust condensation of 30 Mq. The CORE source is resolved into two compact and unresolved sources which 
can be well-fit by two homogeneous hypercompact H II regions each one photo-ionized by a BO. 5 ZAMS star, 
or by free-free radiation from shock-ionized gas resulting from the interaction of a jet/outflow system with the 
surrounding environment. The spatial distribution and kinematics of water masers close to the CORE-N and S 
sources, together with excess emission at 4.5 fim and the detected dust emission, suggest that the CORE source 
is a massive protostar driving a jet/outflow. 

Subject headings: stars: formation — ISM: individual objects (G75.78h-0.34) — ISM: HII regions — ISM: 
dust — radio continuum: ISM 
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1. INTRODUCTION 

Massive stars are key to understanding many physical phe- 
nomena in the Galaxy, however their first stages of forma- 
tion are still poorly understood. One of the main reasons is 
that massive stars evolve more quickly to the main-sequence 
than low mass stars, and they radiate large amounts of ultra- 
violet (UV) photons and drive strong winds while they are 
still deeply embedded and accreting matter (e. g., Beuther & 
Shepherd 2005; Keto 2007). The interaction of the UV radi- 
ation and winds with the surrounding environment, resulting 
in bright radio continuum sources, must be well understood to 
comprehend the formation of high-mass stars. 

Thermal (free-free) radio emission at centimeter wave- 
lengths in regions of massive star formation can have distinct 
origins: i) H II regions photoionized by embedded massive 
stars, with homogeneous density distributions (e.g., Mezger 
& Henderson 1967) or with density gradients (e.g., Olnon 
1975; Panagia & FelH 1975; Franco et al. 2000); ii) clumps 
of gas or circumstellar disks externally ionized by luminous 
early-type stars (e.g., Garay 1987; O'Defl & Wong 1996; 
Zapata et al. 2004); Hi) shock waves arising either in dense 
interstellar gas (e.g., Ghavamian & Hartigan 1998; Araya 
et al. 2009) or from the collision of thermal radio jets from 
young stellar objects (YSOs) with their surroundings (e.g., 
Anglada 1996; Eisloff^el et al. 2000; Rodriguez et al. 2005); 
iv) ionized accretion flows where the material becomes ion- 
ized while accreting onto the massive protostar (e.g., Keto 
2002, 2003, 2007); v) photoevaporated disks where the ra- 
diation of the newly-formed star ionizes and evaporates the 



surrounding disk (e.g., HoUenbach et al. 1994; Lugo et al. 
2004; Avalos & Lizano 2012); vi) equatorial winds with the 
emission produced by small-scale ionized stellar winds (e. g., 
Hoare 2006). All these mechanisms produce radio contin- 
uum sources with spectral indices, cf (5,, oc v"), between -0.1 
and +2 (i.e., thermal emission; see Rodriguez et al. 1989). 
In addition, several works have found negative spectral in- 
dices (typical of non-thermal emission) associated with mas- 
sive YSOs (e.g., Rodriguez et al. 1989; Zapata et al. 2006). 
These non-thermal sources can be young stars with active 
magnetospheres producing gyro-synchrotron emission (e.g., 
Feigelson & Montmerle 1999); synchrotron emission from 
fast shocks in disks or jets (e. g., Reid et al. 1995; Shepherd & 
Kurtz 1999; Shchekinov & Sobolev 2004); or extremely em- 
bedded YSOs where the UV photons from the massive proto- 
stars are heavily absorbed by large amounts of dust, with mass 
column densities >10-' g cm"- (Rodriguez et al. 1993). A 
similar description of these mechanisms producing centimeter 
continuum emission can be found in Rodriguez et al. (2012). 

It is possible — and indeed, probable — that several of 
these emission mechanisms, either thermal or non-thermal, 
occur within massive star formation regions, either simulta- 
neously or at different evolutionary epochs. Any complete 
model of high-mass star formation must address the presence 
of these multiple modes of radio continuum emission. 

The ON-2 star forming complex contains several early-type 
(O and B) stars within a massive (> 10"* Mq, Matthews et al. 
1986; Dent et al. 1988) molecular cloud that spans ~10' on 
the sky. Seen in CO, the cloud has two distinct condensa- 
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Main continuum observational parameters of G75.78+0.34 
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TABLE 2 

Main observational parameters of the VLA and OVRO spectral line observations 
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tions, with a roughly north-south orientation (Matthews et al. 
1986). Multiple ionized regions within the southern conden- 
sation were identified by Matthews et al. (1973, 1977). They 
reported an extended H II region G75.77+0.34, and also the 
ultracompact (UC) H II region G75.78-i-0.34 (hereafter G75); 
the latter is associated with strong OH maser emission first 
reported by Ellder et al. (1969). 

ON-2 lies toward Cygnus-X, and hence it has been prob- 
lematic to assign a reliable distance. Early estimates ranged 
from 0.9 to 5.5 kpc with the nearer distance coming from ex- 
tinction or luminosity arguments and the farther distance com- 
ing from Galactic rotation models and radio recombination 
line velocities. A helpful summary is presented by Camp- 
bell et al. (1982). More recent works, focusing on the UC 
H II region rather than the molecular cloud, tend to adopt 
a kinematic distance of 5.6 kpc (e.g.. Wood & Churchwell 
1989; Hanson et al. 2002). More recently, Ando et al. (201 1) 
observed the water masers associated with the G75 UC H II 
region as part of the VERA (VLBI Exploration of Radio As- 
trometry) project. They measured the trigonometric parallax 
and report a heliocentric distance of 3.83+0.13 kpc. We con- 
sider this distance measurement to be the most reliable and 
we adopt 3.83 kpc for our analysis. This distance places 
G75 close to other star forming regions (G75.76+0.35 and 
AFGL2591; Rygl et al. 2012), and close to the solar circle 
(Ando et al. 2011; Nagayamaet al. 2012). 

At centimeter wavelengths, G75 is dominated by a 
cometary UCH II region reported by Wood & Churchwell 
(1989). Hofner & Churchwell (1996) detected a cluster of 
water masers located about 2" (10"^ AU) southwest of the 
UCH II region, coincident with a compact radio continuum 
source (Carral et al. 1997) with a spectral index of 4-1 .5 + 0.4 
from 6 cm through 3 mm. Franco et al. (2000) modeled this 



compact continuum source as a hypercompact (HC) H II re- 
gion with iig oc r"^. They note that this very steep density 
gradient is probably unrealistic, and mentioned several pos- 
sible causes, including a contribution from warm dust to the 
flux density at 0.7 cm. Additionally, emission from a myr- 
iad of molecular line transitions has been reported in single- 
dish surveys (e. g., Shirley et al. 2003, Roberts & Millar 2007, 
Klaassen & Wilson 2007, Bisschop et al. 2007). Higher angu- 
lar resolution observations of different dense gas tracers show 
that most of the molecular emission comes from the compact 
radio continuum source associated with the cluster of water 
masers (Codella et al. 2010; Sanchez-Monge 2011). Finally, 
up to four distinct outflows have been identified in the ON-2 
cloud core (Shepherd et al. 1997). All this makes G75 an ex- 
cellent target to study the nature of the centimeter continuum 
sources in a massive star forming region. 

In this paper, we present high angular resolution centime- 
ter continuum observations together with 22 GHz H2O and 
44 GHz CH3OH maser observations toward G75. We com- 
plement this data with millimeter continuum and radio re- 
combination line observations, with the goal of deciphering 
the nature of the centimeter continuum sources. 

2. OBSERVATIONS 

2.1. VLA radio continuum observations 

G75. 78-1-0. 34 was observed with the Very Large Array 
(VLAQ) at 6.0, 3.6, 2.0, 1.3, and 0.7 cm from January 1996 
to April 2001, using the array in the CnB, B, and A config- 
urations. In Table [T| we summarize these observations. The 
data reduction followed standard procedures for calibration of 

' The Very Large Array (VLA) is operated by the National Radio As- 
tronomy Observatory (NRAO), a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, Inc. 
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high frequency data, using the NRAO package AIPS. Initial 
images were produced with a robust parameter (Briggs 1995) 
of 1 (see Table [T] for synthesized beams and rms noise levels 
of these images). These data can be grouped into two cate- 
gories based on angular resolution: those with flbeam ^ I'-'O, 
and those with (^eam ~0'.'2. The continuum images at 1.3 cm 
and 0.7 cm, from projects AK490 and AK500, were cross- 
calibrated with the strongest H2O and CH3OH maser compo- 
nents, respectively, which were obser ved s imultaneously with 
the continuum emission (see Section IZ2] and Table |2]l. After 
comparing the initial maps for consistency, we combined the 
Mv-data at the same frequencies to obtain final images with 
better t/v-coverage and sensitivity. The resulting synthesized 
beams and rms noise levels of the combined images are given 
in Table [3] 

At 6 cm, the H II region G75.77+0.34, approximately 1' 
to the southwest (e.g., Riffel & Liidke 2010), produced non- 
imageable artifacts. The shortest baselines (< 5 k/l) were ex- 
cluded, producing essentially no change in the measured flux 
density but significantly improving the quality of the map. 

2.2. VLA HoO and CH^OH maser observations 

The water maser line at 22.23508 GHz (616 -523 transition) 
was observed with the VLA in the A configuration (project 
AK490) simultaneously with the 1.3 cm continuum emission. 
Two different correlator configurations were used, providing 
spectral resolutions of 0.3 and 1.3 km s"'; with velocity cov- 
erages of 21 and 83 km s"', respectively. 

The class I methanol maser line at 44.06941 GHz (7o - 61 
A^ transition) was observed with the VLA in the B configura- 
tion (project AK500) simultaneously with the 0.7 cm contin- 
uum emission. Details of the spectrometer configuration are 
given in Table |2] 

The H2O and CH3OH maser data were calibrated following 
the AIPS guidelines for calibration of high frequency data. 
Self-calibration was performed on the strongest maser com- 
ponent, and the solutions were applied to the spectral line and 
continuum data. The images were constructed using uniform 
and natural weighting to measure the maser positions at the 
highest angular resolution and to estimate the intensity of the 
different maser components, respectively. 

2.3. OVRO observations 

The Owens Valley Radio Observatory (OVRCQ) observa- 
tions at 3 and 1 mm were made in the L (Low) and H (High) 
resolution configurations during September, October and De- 
cember 1997. In March 1998, additional 3 mm observations 
were made in the uH (ultra-High) resolution configuration. 
All the observations were made in the double sideband mode, 
simultaneously observing at 3 and 1 mm. The continuum was 
observed in two 1 GHz channels, centered at 95.78 GHz and 
98.78 GHz for 3 mm and 228.85 GHz and 231.85 GHz for 
1 mm. 

In addition, spectral line modules covered the H40q' 
(99.02296 GHz) and H30q' (231.90094 GHz) radio recom- 
bination lines (RRLs). Each spectral line setup consisted of 
62 Hanning smoothed channels of 2 MHz each, providing 
resolutions of 6.1 and 2.6 km s"' at 3 mm and 1 mm, re- 
spectively. The assumed LSR velocity for both lines was 
-8 km s ' . Bandpass calibration was performed by observing 
the quasar 3C454.3. Amplitude and phase calibration were 

- OVRO was operated by the California Institute of Teclinology witli sup- 
port from the National Science Foundation. 
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Fig. 1. — G75.78+0.34 star-forming region, (a): VLA 6.0 cm continuum 
image. Levels are -3, 3, 6, 10, 15, 25, 35, 45, and 55 times 490 //Jy beam"'. 
Ten-point red stars indicate CH3 OH masers, and three-point yellow stars in- 
dicate H2O masers. The dashed box shows the region zoomed in the bottom 
panel, (b): VLA 3.6 cm continuum image. Levels are -4, 4, 8, 12, 20, 30, 
and 40 times 50 /ijy beam" ' . Synthesized beams of the continuum images are 
shown in the bottom-left comer; spatial scales are indicated in the top-right 
comer of each panel. 



achieved by monitoring BL Lac during the different observ- 
ing tracks. The absolute flux density scale was determined 
from Uranus, with an estimated uncertainty of 20% at 3 mm 
and 30% at 1 mm. The data were reduced using a combi- 
nation of routines from the OVRO reduction package MMA, 
from MIRIAD, and from AIPS. Imaging was performed with 
the task IMAGR of AIPS. In Tablesn]and[2] we list details of 
the continuum and radio recombination line observations. 

2.4. SMA observations 

G75.78+0.34 was observed with the Submillimeter Array 
(SM/fl) in the 1.3 mm (230 GHz) band using the com- 
pact configuration on 2010 June 10. A total bandwidth of 
2x4 GHz was used, covering the frequency ranges 218.2- 
222.2 GHz and 230.2-234.3 GHz, with a spectral resolution 
of ~1 km s"'. System temperatures ranged between 150 and 
250 K. The zenith opacities at 225 GHz were around O.IO 

' The SMA is a joint project between the Smithsonian Astrophysical Ob- 
servatory and the Academia Sinica Institute of Astronomy and Astrophysics, 
and is funded by the Smithsonian Institution and the Academia Sinica. 
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Fig. 2. — Close-up continuum images of the UCHII, EAST and CORE sources. In each panel we show the VLA image at the wavelength indicated in the 
top-right comers. Left column: UCHII source. For all panels, levels are -4, 4, 8, 12, 20, 30, and 40 times the rms noise level of the map: 485, 61, 138, 67, and 
159 /iJy beam"', for 6.0, 3.6, 2.0, 1.3, and 0.7 cm maps, respectively. Center column: EAST source. For all panels, levels are -4, 4, 8, 12, 16, 20, 25, and 30 
times the nns noise level of the map: 61, 138, 67, and 159 fiiy beam"', for 3.6, 2.0, 1.3, and 0.7 cm maps, respectively. Right column: CORE (N and S) source. 
For all panels, levels are -3, 3, 6, 12, 18, 24, and 30 times the nns noise level of the map: 61, 144, 67, and 186 ,uJy beam"', for 3.6, 2.0, 1.3, and 0.7 cm maps, 
respectively. See Table|3]for details of the synthesized beams, and the flux densities and sizes of the sources. Blue crosses in all panels indicate the position of 
the four continuum sources: UCHII, EAST, CORE-N and CORE-S (see Table[3). The scale in arcseconds is shown in the 2.0 cm panels. At the adopted distance 
of 3.8 kpc, 2" con-esponds to 7600 AU while 0.5" is 1900 AU. 



and 0.15 during the 3-hour track. The FWHM of the pri- 
mary beam at 1.3 mm was ~56". Bandpass calibration was 
performed by observing the quasar 3C454.3. Amplitude and 
phase calibrations were made by monitoring 2025-1-337 and 
2015-1-371, with an rms phase of ~40°. The absolute flux 
density scale was determined from Callisto with an estimated 
uncertainty around 15%. Data were calibrated and imaged 
with the MIRIAD software package. The continuum was con- 



structed in the (u,v) domain from the line-free channels. Imag- 
ing was performed using natural weighting, resulting in a syn- 
thesized beam of 6'.'4 x 2'.'7 with a P.A.=-76°, and 1 cr rms of 
8.4 mJy beam"' for the continuum. The two 4 GHz passbands 
include several molecular transitions, including CO, CH3CN, 
and CH3CCH. These molecular line data will be presented in 
a forthcoming paper together with VLA ammonia observa- 
tions (Sanchez-Monge et al., in prep.). 
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TABLE 3 

MULTIWAVELENGTH RESULTS FOR THE YSOs m THE STAR-FORMING REGION G75.78+0.34 
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" Primary beam corrected. Error in intensity is 2a. Error in flux density has been calculated as |^(2(T (Ssourcc/^^bcam)''^ ) + (2crfl„x-scaie) ] , 
where tr is the rms noise level of the map, ^source and ^bcam are the size of the source and the beam respectively, and CTflux-scaic is the error 
in the flux scale, which takes into account the uncertainty on the calibration applied to the flux density of the source (S,, x % uncertainty)- 

'' Coordinates calculated as an average of the coordinates at each wavelength where the source is well detected (i. e., not filtered out). 
Typical discrepancies in the coordinates at different wavelengths are < (/.'l. 



3. RESULTS 

3.1. Centimeter continuum emission 

We detected centimeter radio continuum emission at all 
wavelengths. In Figure [T] we show a global overview of the 
G75 region at 6.0 cm (panel a) and 3.6 cm (panel b). As 
seen in the figure, the centimeter continuum emission is dom- 
inated by three components that we name UCHII, EAST and 
CORE. In Figure |2] we show a close-up of these three sources 
for the combined maps at each frequency, and in Table |3] we 
list their flux densities and sizes, together with the beams, and 



rms noise levels of the combined images. 

The strongest source (UCHII; G75. 7826 -hO. 3429) is a 
cometary UCH II region with an integrated flux density of 
~35 mJy and an angular diameter of ~1" (~0.02 pc at a dis- 
tance of 3.8 kpc; previously imaged by Wood & Churchwell 
1989 and Carral et al. 1997). Located ~6" to the east, we 
identify a compact source (EAST; G75.7830+0.3416), with 
a flux density of ~4 mJy and an angular diameter of ~(y.'2 
(~760 AU at a distance of 3.8 kpc). This source appears 
at the ~5cr level in the 7 mm maps of Carral et al. (1997). 
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Fig. 3. — OVRO and SMA millimeter continuum maps. Levels are -4, 4, 8, 12, 20, 30, and 40 times therms noise level of the map: 0.6, 4, and 8.4 mjy beam"' , 
for OVRO 3.1 mm, OVRO 1.3 mm, and SMA 1.3 mm maps, respectively. See Table[3]and Section |X21 for details of the flux and beam of each image. Blue 
crosses show the position of the radiocontinuum sources, as in Figure|2] 




20 21 ""44^2 44'l 44^0 43^9 

Right Ascension (J2000) 



Fig. 4. — H2O maser distribution in G75. Black contours: VLA 1.3 cm continuum image as in Figure|2]- middle column. Color circles show the position of the 
35 water maser spots (see Table|4). Different colors ai'e used to indicate the maser LSR velocities, according to the color scale (in km s~') on the right hand side 
of the plot. 



Finally, at the head of the cometary arc we find a compact 
source (CORE; G75.7821 +0.3418) slightly elongated in the 
northeast-southwest direction, with a flux density of a few 
mJy and increasing with frequency. This source is coinci- 
dent with the 7 mm continuum source reported by Carral 
et al. (1997) and with the clump of H2O masers (Hofner 
& Churchwell 1996). Our higher angular resolution obser- 
vations (<0'.'3) allow us to resolve the CORE source into 
two distinct (marginally resolved) condensations: CORE- 
N (G75.7821+0.3428) and CORE-S (G75.7820+0.3429; see 
Figure |2]i. 

3.2. Millimeter continuum emission 

In Figure |3] we show the OVRO and SMA millimeter con- 
tinuum images of the G75 region. At 3.1 mm we detect emis- 
sion from all three sources, with the emission from the UCHII 
and CORE sources barely resolved. At 1 .3 mm we only detect 
emission associated with the CORE source, probably due to 
a lack of sensitivity that precludes the detection of faint emis- 
sion from the EAST and UCHII sources. The angular resolu- 
tion of our millimeter continuum images (10-15 times poorer 
than the angular resolution at centimeter wavelengths) does 
not allow us to resolve the components CORE-N and CORE- 
S. In Table[3]we list the flux densities and sizes, and the beams 
and rms noise levels of the OVRO observations. For the SMA 



1.3 mm source, we fit a two-dimensional Gaussian obtaining 
a primary beam corrected flux density of 0.60 + 0.17 Jy with 
a deconvolved size of (4 x 2 ± 1)" at a RA.=(75 ± 10)°. The 
SMA 1.3 mm source, although centered on the CORE posi- 
tion, encompasses the UCHII source and shows an extension 
to the EAST source. In addition, the flux density measured 
in the SMA image is almost twice the flux density measured 
in the OVRO image, suggesting that there is some faint, ex- 
tended emission not detected in the OVRO map. We refrain 
from combining the OVRO and SMA 1.3 mm continuum im- 
ages because of the poorer Mv-coverage and sensitivity of the 
latter. Higher sensitivity maps are needed to confirm the 1 mm 
continuum emission associated with the UCHII and EAST 
sources, and higher angular resolution is needed to identify 
the millimeter emission of CORE-N and S sources. 

At 3.5 mm Riffel & Liidke (2010) measured a flux density 
of 119 mJy with a synthesized beam of 18", while Shepherd 
et al. (1997) measured 75 mJy with a resolution of 5", both 
observations made with the BIMA telescope. Our measure- 
ment (~75 mJy, including all the flux at 3.1 mm) is in good 
agreement with that of Shepherd et al. but somewhat lower 
than that of Riff'el & Liidke. This difference may arise be- 
cause the latter work (with 18 "resolution) was more sensitive 
to weak, extended emission. 
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TABLE 4 

22 GHz WATER AND 44 GHz METHANOL MASERS IN G75.78+0.34 
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Labels H and L (for water masers) refer to the VLA observations 
done with a spectral resolution of 0.33 km s"' and 1.3 km s"', respec- 
tively. Numbers 1^ (for methanol masers) refer to the four groups 
of methanol masers detected (see Figure [T]panel a) 



3.3. HiO and CHiOH maser emission 

Hofner & Churchwell (1996) reported a cluster of water 
masers located 2" southwest from the UCH II region, at the 
same position as the CORE source (as reported by Carral et 
al. 1997). Our H2O maser observations at 22.235 GHz have 
a double aim: to cross-calibrate the radio continuum data at 
1.3 cm, and to observe the masers with higher angular reso- 
lution (-(X'l versus ~0'.'4 of previous observations). Two dif- 
ferent spectral resolutions (~0.3 and ~1.3 km s"') were used 



to observe the H2O maser emission, allowing us to look for 
maser components in a velocity range of (-40, -1-40) km s"'. 
In Figure |4] we show the different H2O maser positions (col- 
ored circles) overlaid on the 1.3 cm continuum image. We 
detected a total of 35 spots with velocities ranging from -36 
to -1-14 km s"'. In Table|4]we list the position, intensity, ve- 
locity, and integrated intensity of each H2O maser component, 
indicating in the last column at which spectral resolution the 
component was detected. The water masers appear to form 
an arc at a distance of ~2" from the head of the cometary 
UCH II region (cf. Figure[T}, with only a few of them directly 
associated with the CORE-S continuum source. We note that 
there are no instrumental offsets to be considered between the 
1 .3 cm continuum and the H2O maser images, since the obser- 
vations were simultaneous. This arc distribution of the maser 
spots was also reported by Ando et al. (2011) with VERA 
observations at 10 mas scales. In Figured we use different 
colors as an indicator of the velocity of the maser component. 
No clear velocity gradients are found in the H2O maser emis- 
sion, although most of them appear slightly redshifted with 
respect to the cloud velocity, wlsr ~ km s~' ( Olmi & Cesa- 
roni 1999; Codella et al. 2010). In Section|52]we discuss the 
possible association between water masers and the continuum 
emission. 

Regarding the class I methanol masers at 44.069 GHz, we 
detect four general locations of emission (with a total of 8 
different components; see Figure [T] panel a). In Table |4] we 
list the position, intensity, velocity, and integrated intensity 
for all the methanol maser components. None of the 44 GHz 
CH3OH spots appear to be directly associated with any of the 
radio continuum sources, but rather are located at a distance 
of 10"-20" (0.2-0.4 pc for a distance of 3.8 kpc) to the north- 
east of the three main YSOs. As in other star-forming regions 
(e. g., Kurtz et al. 2004), class I methanol masers rarely coin- 
cide with other signposts of star formation (e. g., H II regions, 
OH masers, class II methanol masers). Theoretical models 
of methanol masers suggest that the class I masers arise in 
an environment where collisional processes dominate (Cragg 
et al. 1992; Pratap et al. 2008). The G75 star-forming re- 
gion contains up to four different molecular outflows (Shep- 
herd et al. 1997), suggesting that the 44 GHz CH3OH masers 
could be pumped by collisions resulting from the molecular 
outflows. In particular, some of the methanol masers we de- 
tect are close to or within the red lobe of a molecular outflow 
seen in the CO (2-1) and '''CO (2-1) Unes (Sanchez-Monge 
201 1; see also Figure|7]i. Furthermore, the velocities of these 
maser components are red-shifted (see Table |4) which would 
be expected if they are excited by collisions in the red-shifted 
lobe of the molecular outflow. 

3.4. Radio recombination line emission 

We used the OVRO interferometer to observe the radio re- 
combination lines at 3 mm (H40q') and 1 mm (H30q') towards 
G75, to determine the contribution of the ionized gas compo- 
nent at millimeter wavelengths. However, we did not clearly 
detect either of these lines. The spectra are dominated by con- 
tinuum emission coming from either dust or ionized gas (see 
SectionlHi. Our Icr rms noise levels for the line (continuum- 
free emission) are 10 and 15 mJy beam"' for the H40q' and 
H30q', respectively. Although we did not detect the radio re- 
combination lines, the upper limits can be used to constrain 
the emission from ionized gas in the millimeter range (see 
Section l572l for more details). 
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Fig. 5. — Spectral energy distributions for UCHII (top) and EAST (bottom) 
sources. Circles and upper limits correspond to the observational data from 
Table|3] Red dashed lines: linear fit (Sv oc v") to the centimeter data (from 6 
up to 0.7 cm; see Table|5). 



4. FLUX DENSITY DISTRIBUTIONS 

In Figures |5] and |6] we show the flux density distributions 
(FDDs) of the radio continuum sources found towards G75. 
The UCHII source has a flat distribution, with a spectral in- 
dex (a; Sy oc v") of -0.19 + 0.06, typical of optically thin 
free-free emission. This FDD is well-fit by an optically thin 
H II region with an electron density of 3.7 x lO"* cm a di- 
ameter of ~0.019 pc (~3800 AU; consistent with the observed 
deconvolved size), and ionized by a BO type star In Table |5] 
we list the main physical parameters of the cometary UCH II 
region. The emission detected at 3.1 mm and the upper limit 
at 1.3 mm are consistent with the millimeter continuum emis- 
sion coming from ionized gas, suggesting that almost no dust 
emission is associated with the cometary UCH II region, in 
agreement with Riff'el & Liidke (2010). 

The eastern continuum source (EAST), barely resolved 
(cf. Figure |2] and Table O, has a flat spectral index [a - 
-0.16 + 0.12). Its FDD (see Figure|5]l can be fit, at centime- 
ter wavelengths, by an optically thin H II region with diame- 
ter ~0.004 pc (900 AU; consistent with the deconvolved size 
listed in Table [3]) and an electron density of 1.3 x 10^ cm""*, 
requiring the equivalent of at least one BO. 5 star to provide 
the ionizing photon flux. Its physical parameters are listed in 
Table |5] At 3.1 mm there is an excess of continuum emis- 
sion with respect to the optically thin H II region assumption, 
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Fig. 6. — Flux density distributions for CORE (top), CORE-N (middle) and 
CORE-S (bottom) sources. Circles and upper Hmits correspond to the data 
from Table[3] Red dashed lines: linear fit (5v oc v") to the centimeter data 
(from 6 to 0.7 cm; see Table|5]. Blue dotted lines: homogeneous H II region 
with an electron density specified in the panel. Green dotted lines: modified 
black body law for the dust envelope with a dust emissivity index of /J = 1 , 
a source radius of 3", a dust temperature of 50 K, a dust mass of 30 Mq, 
and a dust mass opacity coefficient of 0.9 cm- g"' at 1.3 mm (Ossenkopf & 
Henning 1994). 



probably coming from a dust envelope of ~ 15-36 M© (see last 
column in Table|5]). We refrain from fitting a dust envelope to 
the millimeter points of the FDD, because we have only an 
upper limit at 1 .3 mm. Instead, we estimate an upper limit for 
the dust emissivity index fi < 1 .4, from the 3 mm detection 
and the 1 mm upper limit. 
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TABLE 5 

Physical parameters of the H II regions and dust properties for the sources in G75.78+0.34 







H II region physical parameters" 
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CORE-N 
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> 1.2 X lO' 


> 1.2 X lO'' 


< 1.4 X 10' 


-5 


4.2 X lO''^ 


B0.5 


<40 
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< 30 


CORE-S 


< 200 


> 1.4 X 10' 


> 1.3 X 10'' 
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-5 
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Physical pai'ameters of the H II regions assuming homogeneous ionized gas. For the EAST and UCH II sources the parameters were obtained by fitting 
an optically thin H II region at all centimeter wavelengths, and for the CORE, CORE-N and CORE-S sources the parameters correspond to the fit of an 
homogeneous H II region at all centimeter wavelengths (see Figure|5]. The fit of an homonegeous H II region is obtained from the expression 



2hv^ 



1 



el'y/kn - 1 



(1 - e-^"'">) n« 



(1) 



where Sy is the flux density at frequency v, By(Tc) is the Planck function corresponding to the electronic temperature assumed to be 10 K, nsource 
is the source sohd angle, and r//(y) is the optical depth defined by r//(v) a 0.08235 [ ^^^f^ J [^] ' [j^] ' (Altenhoff et al. 1960). The free 
pai'ameters in the fit ai'e the emission measure and the size of the source. The spectral type is determined from Panagia (1973) using the number of 
ionized photons, N[ . 

S y3,'jj„, and 5 correspond to the (dust) millimeter continuum flux density after subtracting the contribution from free-free emission (from the fits 
shown in Figurejs). For the CORE source, the first value assumes the centimeter emission is described by an homogeneous H II region, and the second 
value assumes the centimeter emission is fitted with a power law fit (S oc v") with a = 1.1. 

Dust and gas mass estimated from the millimeter emission (after subtracting the contribution of the ionized gas). For the UCH II source, all the milhmeter 
continuum emission is thermal ionized gas emission. For the EAST source we assumed a dust emissivity index of 1.5, a dust mass opacity coefiicient of 
0.9 cm- g"' at 1.3 mm (Ossenkopf & Henning 1994), and a dust temperature of 50 and 20 K, respectively. For CORE, CORE-N and CORE-S sources 
we provide the mass estimated assuming a dust emissivity index of 1, a dust mass opacity coefficient of 0.9 cm^ g"' at 1.3 mm, and a dust temperature of 
50 K. The upper limits ai'e due to the low angulai' resolution at millimeter wavelengths that does not allow to resolve the CORE-N and CORE-S sources. 



Finally, we construct the FDD for the CORE source, first 
including all the emission and then for CORE-N and CORE- 
S separately (Figure |6). At centimeter wavelengths, the 
emission is partially optically thick, with a spectral index of 
-1-1.1 + 0.2 for the total emission of the CORE. The major dif- 
ference between the FDDs of CORE, CORE-N and CORE-S 
is the flux density at 7 mm. However, in all cases the centime- 
ter emission can be well-fit by an H II region with an electron 
density of ~10^ cm"-' and a size < 0.001 pc (< 200 AU; see 
Table |5j. An H II region with these properties has a spectral 
turnover in the centimeter regime. At millimeter wavelengths, 
where we spatially resolve the contributions of CORE-N and 
CORE-S, the emission is dominated by hot dust; we estimate 
an envelope mass of 30 Mq and a temperature of 50 K (see 
Table|5). 

5. NATURE OF THE IONIZED GAS EMISSION 

We report the detection of four distinct radio contin- 
uum sources toward the high-mass star forming region G75: 
UCHII, EAST, CORE-N and CORE-S (the latter two cor- 
respond to the CORE source for angular resolutions coarse 
than (X'S). In this section we discuss the nature of these four 
sources. 

5.1. UCHII and EAST: two ultracompact H 11 regions 

The UCHII and EAST sources can be well-characterized as 
homogeneous H II regions with sizes (3800 AU and 900 AU) 
and electron densities (3.7 x lO"* cm"-' and 1.3 x 10^ cm"-') 
characteristic of small ultracompact H II regions (e. g., Kurtz 
2005). They are probably ionized by BO and B0.5 ZAMS 
stars, respectively, and are optically thin at centimeter wave- 
lengths (i.e., spectral index of -0.1; see Table |5]l. Despite 



these similarities, the two sources are distinctly different in 
their infrared emission. In particular, there is substantial IR 
dust emission associated with the EAST source but not with 
the UCHII source. In Figure [T] we show infrared and sub- 
millimeter images of the G75 star-forming complex. Infrared 
emission between 2.2 and 8.0 //m is clearly associated with 
the EAST source, while no infrared emission at these wave- 
lengths is associated with the cometary UCHII source (cf. 
close-up views in Figure |7j. The absence of infrared emis- 
sion from the UCHII source is unusual; ultracompact H II 
regions typically present substantial IR emission (Hoare et al. 
2007). A possible explanation is that the UCH II region may 
be partially obscured by a cloud of gas and dust that absorbs 
the near- and mid-infrared emission. In fact, such a cloud is 
seen in emission at submillimeter wavelengths (cf. emission at 
450 and 850 fim in Figure|7j Di Francesco et al. 2008), and in 
absorption at infrared wavelengths (cf. obscured region in the 
infrared maps in Figure |7), and its peak falls very close to the 
UCHII source. In addition, if the EAST source were located 
closer to the border of the dust cloud (seen in the 450 j^m and 
850 i-im images), its IR emission would be less obscured than 
that of the UCHII source, resulting in the different infrared 
properties of the two sources. The non-detection of UCHII at 
1.3 mm (see Fig. [ 3]l does not contradict this scenario; as we 
noted in Sect. 13.21 the non-dection is probably due to inade- 
quate sensitivity. 

5.2. CORE: hypercompact H II regions or shocks? 

Regarding the emission from the CORE (N and S) source, 
we examine several models to explain the origin of the ionized 
gas emission. Franco et al. (2000) modeled the CORE source 
as an H II region with a density distribution He mo- 
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Fig. 7. — Submillimeter (JCMT) and infrared (Spitzer and UKIDSS) images of G75. Top-left: 450 /jm JCMT continuum image (wliite contours, levels start at 
10% increasing in steps of 10% of the peak intensity 99.6 Jy) overlaid on the 5.8 /im IRAC-Spitzer image. Top-right: 850 /im JCMT continuum image (white 
contours, levels start at 10% increasing in steps of 10% of the peak intensity 10.3 Jy). overlaid on the 4.5 /im IRAC-Spitzer image. The white box shows the 
region zoomed in the bottom panels. The large infrared and submillimeter source toward the southwest of the images con'esponds to the H II region G75.77+0.34. 
A dark region is seen toward the northeast in the infrared images, coincident with the strong submillimeter emission, and associated with the G75 star forming 
complex. The strongest sources in the region, including the large H II region G75.77+0.33, appear saturated (black pixels) in the Spitzer images. Bottom: 
Close-up view of the region studied with interferometers showing the 5.8 fim, 4.5 fim, and 2.2 /jm images. The white contours in the lower left panel show the 
450 fjm JCMT continuum emission as in the top panel. The soHd/dashed white contours in the bottom middle panel show the blue/red-shifted CO (2-1) emission 
(Sanchez-Monge 2011; the contours are 10, 30, 70, 120 and 150 Jy beam"' km s"') tracing the molecular outflow with a direction similar to the elongation 
seen at 4.5 /jm. The white/black stars mark the positions of the radio continuum sources: UCHII, EAST and CORE-N and CORE-S (see Table[3]l. No infrared 
emission is associated with the UCHII source, while strong emission is coincident with the EAST source (see main text for discussion). In all panels, the units of 
the color scales are 10' MJy sr . 



tivated by optically thick centimeter continuum emission at 
frequencies up to ~100 GHz. Our new, higher resolution ob- 
servations suggest that the emission becomes optically thin at 
frequencies ~30 GHz. Figure |6] suggests that a homogeneous 
density distribution can account for the emission from the N 
and S sources. Furthermore, the RRL upper limits are also 
consistent with optically thin emission. Assuming optically 
thin emission for line and continuum, and 3cr upper limits for 
H30q' and H40q' (i. e., 30 and 45 mJy beam"', respectively), 
we can calculate the expected continuum emission from ion- 



ized gas with (Rohlfs & Wilson 2004) 



5l 


= 6940 


Au 


-1 




-1. 15 


V 


Sc 


km s"' 




K 




GHz 



(2) 



where 5l is the line flux density, is the continuum flux 
density, Av is the width of the radio recombination line (as- 
sumed to be 30 km s"'; Kurtz 2005), is the electron tem- 
perature assumed to be 10"* K, and v is the frequency of the 
line. We find that the continuum flux density of ionized gas at 
3 and 1 mm should be <30 and <18 mJy, respectively. These 
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values are in agreement with the homogeneous H II region fits 
given in Section |4] for which the expected flux densities at 3 
and 1 mm are ~3 mJy (i.e., <s;18 mJy). Thus, we conclude 
that the ionized gas emission becomes optically thin for fre- 
quencies >30 GHz. This new result, compared to Franco et 
al. (2000), arises because we resolve the CORE source into 
two components. 

From the FDD analysis (see Section |4|i, both sources 
(CORE-N and S) appear to be hypercompact H II regions 
ionized by BO. 5 ZAMS stars, and separated by a distance of 
0'.'36 (1400 AU at a distance of 3.8 kpc). This scenario sug- 
gests a wide, massive binary system (Sana & Evans 2011). 
We searched the literature for similar systems, with two close 
massive radio continuum sources, and found several exam- 
ples: in G3 1.41 +0.31 there is a system of two centimeter 
continuum sources separated by 0'.'19 (1500 AU at 7.9 kpc; 
Araya et al. 2008); in GlO.47+0.03 there are two similar 
sources separated by 0'.'53 (5700 AU at 10.8 kpc; Cesaroni 
et al. 2010); and in DR21(OH) a similar system is separated 
by 0'.'45 (900 AU at 2 kpc; Araya et al. 2009). Thus, it seems 
that it is not unusual to detect double radio continuum sources, 
separated by ~2000 AU, in massive star-forming regions. 

It should be noted that the radio continuum emission from 
some of the previously listed binary systems is not always in- 
terpreted as arising from photo-ionized H II regions. Araya 
et al. (2009) propose that the radio continuum emission in 
DR(21)OH is free-free radiation originating in interstellar 
shocks, following the theoretical development of Ghavamian 
& Hartigan (1998). In this scenario, the free-free emission 
arises from shock-ionized, rather than photo-ionized, gas, be- 
ing comparable to the Herbig-Haro objects seen in the optical. 

We note that the infrared emission at 4.5 yum (see Figure|7]i 
shows an elongated structure 'emanating' from the position 
of the CORE source. The fluxes of this structure in the four 
IRAC/Spitzer bands are 4 mJy at 3.6 //m, 52 mJy at 4.5 fim, 
96 mJy at 5.8 jim, and 115 mJy at 8.0 fim, and thus it fulfills 
the criteri^ defined by Chambers et al. (2009) to be classi- 
fied as a 'green fuzzy' (see also Cyganowski et al. 2008). The 
excess 4.5 fim emission can be produced by shocks associ- 
ated with outflows, scattered continuum in outflow cavities, 
obscuration affecting the emission at 3.6 fim, or fluorescence 
Ht emission (e.g., Noriega-Crespo et al. 2004; Smith et al. 
20"06; Qiu et al. 2008; De Buizer & Vacca 2010; Takami et al. 
2010; VaiTicatt 2011; Simpson et al. 2012; Lee et al. 2012; 
Takami et al. 2012), and only spectroscopic studies in the 
mid-infrared can clearly determine the origin of the 4.5 fim 
excess. However, in G75, the northeast-southwest direction 
of the 4.5 jum elongated structure is consistent with the direc- 
tion of the molecular outflow likely associated with the CORE 
source (see the solid/dashed contours in the bottom middle 
panel of Figure |7] Sanchez-Monge 2011) which is also con- 
sistent with the direction of the outflow reported by Riffel & 
Liidke (2010) at larger scales, thus favoring the interpretation 
of the 4.5 yum emission arising from shocks. In addition, the 
spatial distribution of water masers (Figure |4) is also indica- 
tive of shocks and jets, as found in other star forming regions 
(e. g., IRAS 20126+4104: Hofner et al. 2007). A mifliarcsec- 
ond kinematical study of water masers in G75.78+0.34 (see 
Figure 3 in Ando et al. 201 1) is consistent with a shock inter- 

* Chambers et al. (2009) classified an extended mid-infrared source as a 
'green fuzzy' according to the following color criteiia: (1) the 4.5 fim to 
3.6 fjm flux ratio is > 1 .8, (2) the 4.5 /im to 5.8 pm flux ratio is >0.40, and (3) 
the 4.5 /jm to 8.0 /jm flux ratio is >0.45. 



pretation for the centimeter continuum emission. The group 
of masers associated with CORE-S have velocities mainly 
in the southwest direction (similar to the elongated structure 
seen at 4.5 yum), while the two groups of maser spots located 
to the northwest and southeast of CORE-S show a velocity 
field suggesting expansion. Together with the arc morphol- 
ogy of all the maser spots this is suggestive of a bow-shock, 
implying that at least the CORE-S emission might arise from 
shocks. 

A stellar jet/outflow with a mass-loss rate larger than 
4 10"* Mq yr"' and a terminal velocity of order several hun- 
dred km s"', shocking the ambient medium, would have suf- 
ficient energy to ionize, heat, and move the gas of both the 
CORE-N and CORE-S sources (see Appendix IaIi. Although 
the cooling time (and hence recombination) of the high den- 
sity gas is very short (Franco et al. 1994), the jet/outflow 
would continuously re-ionize the material. In this case, the 
exciting source would be located between the two radio con- 
tinuum sources, with its own H II region quenched by the 
high density medium. I.e., a jet/outflow with such a high 
mass-loss rate would be neutral because the ionization front 
would be trapped very close to the star. The possibility of a 
neutral jet/outflow shocking the ambient gas cannot presently 
be ruled out because our observations are sensitive to ionized 
gas. This scenario can be tested with future centimeter con- 
tinuum observations with the JVLA, to determine if CORE-S 
(and the water maser spots) have been displaced with respect 
to the other continuum sources. Assuming a jet velocity of 
order several hundred km s \ relative motions could be de- 
tected within the next five years. 

As noted in Section [1] other mechanisms might also cause 
the centimeter emission from CORE-N and S. However, we 
consider most of these to be unlikely candidates. Equato- 
rial winds (Hoare 2006) or thermal radio jets (Anglada 1996) 
typically predict constant spectral indices of +0.6 and a rela- 
tion between the source size and the frequency. Our obser- 
vations indicate spectral indices — hi; the uncertainties in the 
deconvolved source sizes are too large to determine a size- 
frequency dependence. A photoevaporated disk model has 
been proposed to explain the continuum emission of double- 
peaked sources with an hour-glass morphology (Lugo et al. 
2004). Although the CORE N and S peaks might be fit as a 
photoevaporated disk, recent studies indicate that the peaks 
from photoevaporated disks should change their separation 
with frequency, moving closer at higher frequencies (Tafoya 
et al. 2004). The peak positions of N and S do not change 
between 3.6 and 0.7 cm to a precision of ~0'.'02. 

Summarizing, the centimeter continuum emission from 
CORE N and S can be well fitted by homogeneous H II re- 
gions, each one photo-ionized by a BO. 5 ZAMS star Alter- 
natively, free-free radiation from shock-ionized gas resulting 
from the interaction of a jet/outflow system with the surround- 
ing environment appears to be a viable scenario as well. Other 
possible mechanisms are deemed unlikely. 

6. CONCLUSIONS 

We have carried out subarcsecond angular resolution obser- 
vations with the VLA in the centimeter continuum and in H2O 
and CH3OH maser emission toward the massive star forming 
complex G75.78+0.34. We have complemented these obser- 
vations with OVRO and SMA millimeter continuum and radio 
recombination line observations, and submillimeter/infrared 
data from on-line databases. Our conclusions can be summa- 
rized as follows: 



12 



Sanchez-Monge et al. 



1 . Our radio continuum data reveal centimeter continuum 
emission at all wavelengths, with the emission com- 
ing from four distinct compact sources: UCHll, EAST, 
CORE-N and CORE-S. The two last sources appear un- 
resolved (CORE source) when observed with angular 
resolutions >0'.'3. 

2. The strongest source, UCHll, is a cometary UCH II re- 
gion with a size of 0.02 pc (3800 AU), an electron den- 
sity of 3.7 X 10"* cm"^ and is ionized by the equivalent of 
a BO ZAMS star The EAST source, located ~6" to the 
east of UCHll, is also an UCH II region with an electron 
density of 1 .3 x 10^ cm"-', ionized by a BO. 5 equivalent 
ZAMS star, but with a smaller (barely resolved) size of 
0.004 pc (900 AU). The millimeter continuum emission 
associated with these two UCH II regions has different 
origins: for the UCHll source it probably traces the ion- 
ized gas, while for the EAST source there is a millime- 
ter excess, suggesting that this source is still embedded 
in a compact dust clump also detectable at mid-infrared 
wavelengths, with a mass of 15-36 M©. We suggest 
that the non-detection of thermal dust emission in the 
UCHll source results from foreground extinction, not 
from an absence of warm dust. 

3. The CORE-N and S sources are located 2" to the south- 
west of the cometary arc of the UCHll source, and 
are associated with water maser and dense molecular 
gas emission. The two sources are very close to one 
another (0'.'36, 1400 AU) and have similar properties: 
flux density increasing with frequency, and unresolved 
emission. The free-free continuum emission may be 
produced by two hypercompact H II regions ionized 
by BO. 5 ZAMS stars or by the collision between a 
jet/outflow and the surrounding environment. A dust 
clump of ~ 30 M0 associated with the CORE source is 
detected at millimeter wavelengths. 

4. We have also reported high angular resolution observa- 
tions of H2O and class 1 CH3OH maser emission. The 
class 1 methanol maser spots are found far from any of 
the radio continuum sources, and are probably excited 



by collisional processes due to the presence of multi- 
ple molecular outflows. The water maser spots appear 
close to the CORE-N and S sources, with only a few 
of them spatially coincident with the S source. The 
spatial distribution of the water masers and the kine- 
matics observed at milliarc second angular resolution, 
suggest that they could originate in a shock between 
a jet/outflow and the local surroundings. In addition, 
extended 4.5 jum emission likely tracing shocked gas 
could be associated with a jet driven by the CORE 
source. 

In summary, we have characterized the radio-continuum 
emission toward the massive star forming complex 
G75.78+0.34, revealing the presence of four sources: 
two ultracompact H II regions (named UCHll and EAST) 
ionized by BO-BO. 5 ZAMS stars, one with an extended and 
cometary shape and the other being barely resolved; and two 
compact sources (named CORE-N and CORE-S) associated 
with at least one massive protostar embedded in a dense and 
massive envelope and with hints of driving an outflow. 
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APPENDIX 



ENERGY BUDGET 



Consider a star that ejects into a cloud a bipolar jet/outflow with a mass-loss rate Mj and a jet velocity Vj. The jet will collide 
with the ambient material in two opposite directions and will produce two shocked ionized regions. The luminosity of the jet, Lj, 
is given by 

i2 



Lj^^Mjvf ^ 7.9 X 10 



35 



Mi 



10-5 Mq yr- 



500 km s" 



erg s 



(Al) 



The shocked ionized gas is post shock ambient gas pushed and heated by the stellar bipolar jet. For simplicity, we will assume 
both shocked ionized regions have the same physical characteristics. Then, the minimum luminosity required to keep the post 
shock gas ionized is 



Li = LoM = 2.2 X lO"* 



1046 S-1 



erg s 



(A2) 



where £0 = 13.6 eV is the ionization energy of the hydrogen atom and A^i is the rate of ionizing photons inferred from the radio 
continuum emission of the ionized region. Furthermore, the energy lost by radiation from each region is 



A{T) Vi = 4.2 X 10 



34 



lO*' cm-3 



100 AU 



erg s 
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where we assume a cooling rate per unit volume A(10'* K)/(«e Wp) ~ 3 x 10 erg cm-' s ' (Osterbrock 1989), a volume of each 
region with radius Vi - 47rrj'/3, and assume that the electron density is equal to the proton density, - «p. Also, we estimate 
the kinetic luminosity of each ionized region averaged over its lifetime as 

Lk ~ = 1.5 X 10^4 

^ Across 

where the mass of ionized gas is M\ - inn Vi, the velocity is i>i, and Taoss is the crossing time that can be estimated from 
Tcioss = Ijvi, where / is the distance from the ionized region to the star emitting the jet. 

We consider now the particular case of the CORE-N and CORE-S sources. From Table|5]we take «i ~ lO*" cm"^, n ~ 100 AU, 
and A^i ~ 4 X 10"*^ s '. We assume that each core is the ionized working surface (WS) of the jet moving into the cloud. For strong 
shocks, the speed of each WS is given by ram pressure balance as v\ - /3vj/(l +/3), where /3 = ^Jpi/p■d is the square root of the ratio 
of the jet and ambient densities (e. g., Raga et al. 1990). For large jet mass-loss rates, close to the emitting source, we assume 
j6 ~ 1. Then, for a jet velocity yj ~ 500 km s"' (of the order of the escape speed of a 10 star assumed to be the source of the 
jet) the shocked ionized gas moves with a speed, Ui ~ u^/l ~ 250 km s"'. Because the two continuum sources are separated by 
~ 1700 AU, we assume a distance to the jet source in between CORE-N and CO RE-S, / ~ 850 AU . Then, the resulting crossing 
time for both ionized regions is Tnoss ~ 16.2 yr Introducing these values in Eq. IA2I IA3I and IA4I we obtain for each shocked 
ionized region Li ^ 8.8 x 10"''' erg s"', L,-ad ~ 4.2 x lO"'"' erg s"', and Lk ~ 1-5 x 10'''* erg s"'. 

Finally, one can now estimate the minimum mass-loss rate and velocity of a jet necessary to ionized, heat, and move the 
CORE-N and CORE-S sources. The energy budget has to fulfill the relation Lj > 2(Li + Lrad + ^k)- From Eq. lAll the minimum 
mass-loss rate of the bipolar jet/outflow that heats and moves the ionized the regions would be Mj > 3.7 x 10"^ Mq yr~'. 
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